import pyproj


# 定义坐标系统
def zhuan(x,y):
    # 定义坐标系统 - 兰勃特投影，中央子午线102°，原点纬度30，割线一30，割线二60
    lon_0 = 102.0
    lat_0 = 30.0
    lat_1 = 30.0
    lat_2 = 60.0

    lb = pyproj.Proj('+proj=lcc +lat_1={} +lat_2={} +lat_0={} +lon_0={}'.format(lat_1, lat_2, lat_0, lon_0))
    # 将兰勃特投影坐标转为WGS84坐标系
    wgs84 = pyproj.Proj("+init=EPSG:4326")
    # 将兰勃特投影坐标转换为WGS84经纬度坐标
    lon, lat = pyproj.transform(lb, wgs84, x, y)
    ncdatadict = dict()
    ncdatadict['lat'] = lat
    ncdatadict['lon'] = lon
    return ncdatadict


if __name__ == '__main__':

   a =  zhuan(2000,1000)
   print(a['lat'])
